Fidelity approach to the disordered quantum XY model 
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We study the random XY spin chain in a transverse field by analyzing the susceptibility of the 
ground state fidelity, numerically evaluated through a standard mapping of the model onto quasi- 
free fermions. It is found that the fidelity susceptibility and its scaling properties provide useful 
information about the phase diagram. In particular it is possible to determine the Ising critical line 
and the Griffiths phase regions, in agreement with previous analytical and numerical results. 



Introduction. - In the last few years concepts borrowed 
from quantum information theory have proven useful in 
characterizing the critical behavior of quantum many- 
body systems In particular, a geometric approach 
to the study of quantum phase transitions (QPTs), i.e. 
the fidelity analysis, has been shown to be an effective 
way of characterizing distinct phases of quantum systems 
@, Ell BUSH- Previously, the fidelity approach has 
been applied to a variety of homogeneous systems. In 
this work we extend these studies to disordered quan- 
tum systems. Specifically, we investigate the behavior 
of the fidelity susceptibility of the disordered XY model 
in a transverse field. It is well known that the presence 
of quenched disorder can have drastic effects on critical 
properties. The appearance of new universality classes 
and novel states of matter such as the Griffiths phase 
are two important examples d, [l(| UH, H2| ■ The aim of 
the present work is to show what can be inferred about 
the physics of the disordered quantum system from the 
properties of the fidelity susceptibility. 

The Hamiltonian of the disordered XY chain is given 

by 
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where a^^'^ 1 are Pauli spin matrices, and 7$ and Xi 
are sets of independent random coupling and field vari- 
ables with distributions 7r(7j) and p(Xi). Note that due 
to gauge symmetry the Hamiltonian |T]) can be chosen 
to have only positive couplings and fields. This model 
can be mapped onto a system of quasi-free fermions with 
periodic boundary conditions, and an exact expression 
for the fidelity susceptibility is obtained which depends 
explicitly on the random parameters characterizing the 
ground state of the system. In this work, we investigate 
the statistical properties of the fidelity susceptibility 0, [E| 
for relevant regions of parameter space. 

Gaussian distributions are used for the random vari- 
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where Xi is either the field or the coupling at position i 
on the chain, x is the respective average value and a 2 is 
the variance. 

Previous Results. - The pure XY chain has been ana- 
lytically solved in In the absence of disorder two 
different quantum phase transitions are present. Follow- 
ing the standard notation, we refer to the QPT driven 
by the transverse magnetic field A as the Ising transi- 
tion, and to the QPT driven by the coupling parameter 
7 as the anisotropy transition. The Ising transition sepa- 
rates a ferromagnetic ordered phase from a paramagnetic 
quantum-disordered phase, whereas the anisotropy criti- 
cal line is the boundary between a ferromagnet ordered 
along the x direction and a ferromagnet ordered along 
the y direction. 

A major improvement in the understanding of the ef- 
fect of disorder on the physics of quantum magnets has 
been achieved with the use of the strong-disorder renor- 
malization group technique (SDRG) by Dasgupta and 
Ma [3, and further developed by Fisher [nf [nj]. The 
correctness of this method has been corroborated both 
by numerics [HI, [Hj] and analytic exact studies [13, EH • 
In the work of McKenzie and Bunder [13, EH the critical 
behavior of the disordered XY chain in a transverse field 
has been studied using a mapping to random-mass Dirac 
equations. The properties of the solutions of these equa- 
tions imply the disappearance of the anisotropy transi- 
tion in the presence of disorder. Furthermore, Griffiths 
phases are predicted to appear both around the Ising 
critical line and the anisotropy 7 = line. 

For 7=1 the XY random chain is closely related to 
the random transverse-field Ising chain (RTFIC), which 
is a prototypical model for disordered quantum systems. 
Since it is representative of the universality class of Ising 
transitions for all values of 7, let us briefly review what 
is known for this model. The Hamiltonian of the RTFIC 
is H = — Y,iJo [Ji a f °f+i + ^°f] ) where J t and h l are 
random couplings and fields respectively. The system is 
at criticality when the average value of the field equals 
the average value of the coupling. Using the SDRG one 
obtains that, at the quantum critical point, the time scale 
r and the length scale L are related by lnr ~ L 1 / 2 . 
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This results in an infinite value for the dynamical expo- 
nent z at criticality. The distribution of the logarithm of 
the energy gap, lne, at criticality broadens with increas- 
ing system size, in accordance with the scaling relation 
lne ~ —L 1 / 2 In the vicinity of the critical point 

the distribution of relaxation times is broad due to Grif- 
fiths singularities. This region of the parameter space, 
the Griffiths phase, is characterized by a dynamical ex- 
ponent z which depends on the distance from the critical 
point. This dependence is one of the hallmarks of the 
Griffiths phase. 

Method. - The main idea of the fidelity approach is to 
detect QPTs through enhanced orthogonalization rates 
between ground states ^(a;)) nearby in parameter space. 
The orthogonalization is signaled by a drop in the fidelity, 
F(x,x + Ax) = |($(a:)|\&(a: + Ax))\, at the critical point. 
The fidelity susceptibility is a related quantity with a 
more transparent physical meaning and whose be- 

havior is more suitable for numerical analysis. It is de- 
fined as 
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In [j| it was shown that x is related to the dynamic struc- 
ture factor of the relevant operator associated with the 
transition. A generalization of this result, valid for the 
so-called geometric tensor, has been given in [4|. 

Previous works have characterized the pure XY spin 
chain using the fidelity approach @, 0, 0] and the 
quantum Chernoff bound [2l| . The mapping of the spin 
model onto the quasi- free fermion Hamiltonian [l3j . 



L 



H = > clAac 



I ±{4 



Bye] + B.C. 



(4) 



yields an explicit BCS-like form for the ground state 
|*) = Wexp (i Ef,k=i c j G ikct) |0), where N is a nor- 
malization factor. 

The fidelity of the ground states evaluated at slightly 
different parameter values (coupling or magnetic field) x 
and x + Ax has a simple analytical expression. Defining 
the matrix Z(x) = A(x) — B(x) and the unitary part of 
the polar decompositions of Z(x) and Z = Z(x + Ax) as 
T(x) and T = T{x + Ax), respectively, the fidelity can 
be written as 
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Note that the matrix G defining the ground state is sim- 
ply the Cayley transform of T [lj| . 

In the following, we will use an alternative expression 
for the fidelity susceptibility, obtained in the limit of 
small Ax, 
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with || ■ \\p the Frobenius norm. Eq. ([6]) is obtained from 
((5|) via standard algebra. We have numerically evaluated 
the fidelity susceptibility using (|6|) for relevant regions 
of parameter space of the disordered XY model. The 
numerical analysis has been performed on two sets of 
system sizes, i.e. {128, 256, 512}, and {400, 410, 500} 
in steps of 10. We have taken 50,000 disorder realizations 
for all sizes except for those larger than 400, in which case 
we used 10,000 realizations. 

Results. - We consider the Hamiltonian |TJ), where the 
couplings 7^ and the transverse fields A,; are independent 
random variables with Gaussian distributions centered 
around A = [Ai] avc and 7 = [7i] avo , both with standard 
deviation a. [-] avc denotes the arithmetic mean over the 
disorder realizations. 

A scaling analysis has been performed using argu- 
ments first developed in [4j|. Following that reference, 
we can express the fidelity susceptibility as an integral 
in imaginary time x — J^ oc> < ^ ttG ( t )j where G(r) = 
0(t) -c d x H(T)d x H(0) ^S> is the connected correlation 
function of the conjugate operator in the Hamiltonian 
associated with the driving parameter in the transition, 
and 8 is the Heaviside step function. For example, in 
the case of the Ising transition we have G(r) = 9{t) <C 
Yli j a i ( T ) a j (®) The average fidelity susceptibility 
can then be written as [x] ave = drr [G(r)] avc - L A * , 
where the finite-size scaling dimension A v of IYL,„ is 
given by A x = 2z + 2 - 2A [4|. A is the scaling di- 
mension of the QPT conjugate operator (V_\ erf in the 
case of the transition driven by A), and in general Ao 
depends on the parameters 7 and A. 

For the XY chain without disorder, in the quantum 
critical regions x scales as x ~ L 2 , whereas away from 
the critical region x ~ L. Since for finite system sizes the 
quantum critical region has a finite width, A x is 1 for all 
but a narrow range of A (or 7), having a maximum of 2 
for A = 1 (7 = 0). 

With this disorder-free behavior in mind, we now study 
Wave about the Ising transition, driven by the coupling 
A. In our numerical studies we have focused on the case of 
the RTFIC, where j = 1. Qualitatively all of our results 
on the critical behavior of the fidelity susceptibility hold 
true for other values of 7, since the universality class of 
the model does not change in the range 7 S (0, 1]. 

Fig. QJa) shows [x] avc as a function of A for the clean 
case and a = 0.1,0.3 disorder strengths. The averaged 
fidelity susceptibility displays a local maximum at the 
Ising critical point, which for the disordered case is 
shifted slightly away from the clean value of A = 1 
due to finite size effects. Fig. [IJb) shows A X (A), the 
finite-size scaling dimension of [x] a vc f° r ^he same 
set of disorder strengths. The disorder leads to a 
broadening in the peak of A x , which is consistent with 
the presence of a Griffiths phase. Note that far from 
the Ising critical point [x] ave scales strictly extensively, 
while in the vicinity of the critical point the scaling 
becomes superextensive. For the weaker noise this 
scaling is nearly quadratic, as in the clean case, while 
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with stronger noise the maximum scaling dimension is 
correspondingly reduced. Qualitatively, the reduction 
of the maximum scaling dimension may be ascribed to 
the presence of rare regions whose extent effectively 
determines the critical behavior. The linear extension 
of rare regions is smaller than the overall system size 
determining the critical behavior in the clean case. 
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FIG. 1: Ising transition at 7 = 1. (a) Average fidelity sus- 
ceptibility, [x] ave , for L — 512 and 10 4 realizations, with 
varying disorder strengths a € {0,0.1,0.3}, (b) the associ- 
ated finite-size scaling dimension A x of [x] avc , ( c ) probability 
distribution of lnx at the Ising transition for system sizes 
L — 128,256,512 and disorder a — 0.1, (d) distribution of 
In x away from the Ising transition for the same disorder and 
range of system sizes. 

In Fig. [TJc) and (d) we plot the distribution of 
the fidelity susceptibility over many realizations at the 
Ising critical point and away from it, for system sizes 
L = 128, 256, and 512. We choose to plot the distri- 
bution of In x instead of x itself because, in analogy to 
other physical quantities, the presence of disorder greatly 
broadens the distribution. As the system size increases, 
note that the probability density function of In \ broad- 
ens for A = 1, but becomes narrower away from critical- 
ity. Indeed, this broadening behavior persists for a range 
of values of A about the critical point. This is typical of 
disordered systems, and is analogous to the absence of 
self-averaging of some physical observables. 

The Griffiths phase around the Ising critical point can 
be detected by looking at the scaling dimension of the 
fidelity susceptibility and at the properties of the dis- 
tribution of lnx, in accordance with the relation A x = 
2z + 2 — 2Ao- The following analysis of the region about 
the anisotropy line further supports this conclusion. 

Although for the disordered XY model the 7 = line is 
not critical, as it is in the pure case [UUII) the presence 
of Griffiths singularities still has highly non-trivial effects 
on the fidelity susceptibility in the vicinity of 7 = 0, as 
shown in Fig. HJa). Specifically, in the presence of disor- 
der the peak in [x] ave (7) splits into two peaks, symmet- 
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FIG. 2: Elimination of anisotropy transition at 7 = due 
to disorder, (a) Average fidelity susceptibility [x] avc i near 
7 = for L = 500, A = 0.2, a = 0.1, and 10 4 realizations, (b) 
Finite-size scaling dimension of [x] avc (in this case 7 = 0.036), 
(c) probability distribution of In x at 7 = and at the value 
of 7 corresponding to the maximum of [x] ave ' (d) distribution 
of In x far away from the anisotropy line, where the finite-size 
scaling is extensive. 
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FIG. 3: Average fidelity susceptibility for various disorder 
strengths a € {0, 0.1, 0.2, 0.3}, with fixed system size L = 400. 
Here A = 0.5, and the derivative in Eq. (((J is taken along 7. 



rical about 7 = 0. Note that 7 = is a special case for 
the XY chain. With zero anisotropy and noise only in 
the field A^, Bunder and McKenzie [18| showed that the 
density of states does not diverge at zero energy. This 
implies that this point is not critical and does not belong 
to a Griffiths phase. We observe similar behavior with 
disorder in both the field and anisotropy. At 7 = the 
fidelity susceptibility scales only extensively, as it does 
in the non-critical region. This suggests that the 7 = 
point does not belong to the Griffiths phase, having char- 
acteristics of a point which is non-critical and away from 
any Griffiths phase. This is further corroborated by the 
non-monotonic dependence on 7 of the associated scal- 
ing dimension shown in Fig. [2jb). At 7 = 0, one finds 
A x = 1, whereas in the interval < 7 < 0.075 the scaling 
dimension A x (7) exhibits a non- universal dependence on 
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the driving parameter 7, indicating the presence of a Grif- 
fiths regime. Note that the observed maximum is not to 
be seen as an indication of a QPT. Rather, it originates 
from the competition between the scaling properties of 
X in the Griffiths phase and at the 7 = line. In Figs. 
EJc) and (d), we show P(ln%) at 7 = 0, at the point 
where x{l) an d ^x(7) both peak, and far away from the 
anisotropy line. In analogy to the Ising transition, the 
probability distribution function in the Griffiths regime 
is broad and asymmetric due to absence of self-averaging, 
whereas far away from it its shape is symmetric and its 
distribution is much narrower. To complete the discus- 
sion of the effects of disorder on the anisotropy transi- 
tion in Fig. [3] we plot the average fidelity susceptibility 
for a fixed system size and various disorder strengths, 
including the clean case. Notice that as the disorder 
strength is increased, the original peak disappears and 
the new maxima in the fidelity susceptibility are sym- 
metrically located around the 7 = line, at a distance 
which increases with disorder. Much like the Ising tran- 
sition, the maximum value of the fidelity susceptibility 
decreases with increased disorder. We believe that this 
can be explained again in terms of the extension of rare 
regions. 

Conclusions. - In this work we have applied the fidelity 
approach to the study of the disordered XY chain in an 
external magnetic field. We have found that the fidelity 



susceptibility is able to provide the phase diagram for 
this model. In the case of the Ising transition, we ob- 
tain results which are consistent with what is already 
known in the literature. In the parameter region around 
the 7 = line the scaling analysis of the fidelity sus- 
ceptibility shows the disappearance of the QPT and the 
emergence of a Griffiths phase, in accordance with similar 
analytical and numerical results. As far as we know, this 
result has not been obtained before for this distribution 
of disorder both in the couplings and in the fields. This 
is nontrivial, since it is known that choosing a different 
parametrization for the disorder can modify the critical 
behavior 
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We plan to further investigate the relevance of disor- 
der on the fidelity susceptibility in future works. Other 
aspects that will be studied with more details are the ex- 
tent of the Griffiths phase together with its dependence 
on disorder strength and the probability distribution of 
disorder. 
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